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ABSTRACT 


This  work  deals  with  the  time-dependent  quadrupole 

7  0 

interaction  of  Ga  in  a  Zn  lattice  for  the  879  keV  level 
with  spin  4. 

2 

The  quadrupole  frequency  has  a  value  w  =  e  q Q / h  = 

315  MHz  and  it  changes  (+3  ±  6)%  from  300  K  to  600  K.  The 
derived  quadrupole  moment  is  0.5  b  when  the  field  gradient 
is  estimated  according  to  the  proposed  universal  correlation 
between  the  local  and  lattice  part  of  the  field  gradient.  The 
temperature  dependence  might  indicate  a  larger  quadrupole 
moment. 


The  measured  anisotropy  coefficients  are  consistent 
with  a  CN-model  for  spin  4  and  the  published  anisotropies 
of  the  decay  of  the  1102  keV  level  with  spin  4~. 

The  ti me- i ntegra ted  values  of  the  attenuation 
coefficients  give  in  combination  with  the  unperturbed  values 
the  indication  that  the  assumption  of  a  static  electric  field 
gradient  is  probably  not  correct  for  all  nuclei  at  higher 
temperatures . 

The  measurements  do  not  yield  a  definite  value  of  the 
spin  of  the  879  keV  level  or  of  the  asymmetry  parameter  n. 
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CHAPTER  I 


INTRODUCTION 

When  a  nucleus  in  an  excited  state  is  formed  in  a 
reaction,  the  population  of  the  magnetic  substates  is  not 
necessarily  equal,  but  the  nucleus  may  be  aligned.  This 
causes  the  decay  of  the  nucleus  to  be  anisotropic.  The 
degree  of  anisotropy  yields  information  about  the  spins  and 
of  the  levels  involved  in  the  decay.  When  the  excited  state 
formed  has  an  electric  or  magnetic  multipole  moment  and 
experiences  an  appropriate  field  for  this  moment,  it  will 
precess.  This  is  the  well-known  Larmor  precession  for  the 
case  of  the  magnetic  dispoles.  This  precession  will  also 
rotate  the  angular  distribution  of  the  decay  of  the  level. 

To  see  the  effect  of  the  precession  the  level  has  to  have 
a  lifetime.  The  technique  which  looks  at  the  time-dependence 
of  this  precession  in  times  longer  than  five  nanoseconds  is 
called  time-dependent  perturbed  angular  distribution  method. 

This  thesis  describes  studies  of  pertubation  of  y  ray 
angular  distributions  due  to  the  interaction  of  the  nuclear 
quadrupole  moment  and  the  electric  field  gradient  (E.F.G.) 
at  the  nuclear  site. 

The  interaction  is  characterized  by  the  spin  of  the 

p 

level  and  the  quadrupole  frequency  e  qQ/h,  where  eq  is  the 
largest  component  of  the  field  gradient  and  Q  is  the  quadrupole 
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moment  of  the  nucleus.  All  three  parameters  are  obtainable 
from  a  time-dependent  angular  correlation  or  distribution 
measurement  (TDPAC  or  TDPAD).  In  the  correlation  experiment 
the  nucleus  is  aligned  by  the  detected  first  decay,  which 
feeds  the  level,  and  the  angular  distribution  of  the  second 
decay,  depopulating  the  level,  is  measured.  The  distribution 
experiments  involves  the  nucleus  aligned  in  a  reaction. 

TDPAC  is  the  most  applied  technique  of  these  two,  for  TDPAD 
experiments  require  a  pulsed  beam  to  define  the  time  of 
creation  of  the  level,  while  TDPAC  uses  the  populating  decay 
as  time  reference  point.  University  of  Alberta  has  a  pulsed 
beam  with  a  width  of  1  ns  and  the  level  of  interest  is  not 
populated  by  a  long-lived  parent,  TDPAD  was  the  way  to  perform 
these  experiments. 

Since  the  frequency  is  a  product  of  two  unknown 
quantities,  i.e.  q  and  Q,  the  results  of  one  measurement  don't 
yield  definite  values  for  both  quantities.  By  looking  at  the 
behaviour  at  various  temperatures  or  in  different  hosts  or 
compare  the  value  of  the  frequency  with  those  of  other  ions 
in  the  same  host,  one  is  able  to  obtain  more  information  about 
the  various  contributions  to>  the  field  gradient,  i.e.  the 
lattice  part,  the  local  part  and  radiation  damage.  The  other 
possibility  to  calculate  one  value  and  deduce  the  other  is 
up  to  the  present,  in  most  cases,  impossible. 

Three  sets  of  measurements  were  carried  out.  The 
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first  concerned  the  interaction  of  Ga  in  a  metallic  Zn 
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target.  In  the  study  of  the  spin  and  parities  of  levels  in 
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the  isotope  Ga  there  appeared  a  discrepancy  between  the 
data  obtained  from  particle  work  and  the  y  ray  anisotropy 
of  the  879  keV  level  ( D  o  7  2 ) .  Since  this  level  has  a  halflife 
of  23  ns,  and  the  Zn  target  used  in  the  y  ray  study  is  non- 
cubic  there  is  the  possibility  of  a  pertubation  of  the 
anisotropy  due  to  a  quadrupole  interaction.  The  other 
collaborators  of  the  published  paper  (Hu75)  had  measured  the 
time  dependence  of  the  anisotropy  at  room  temperature.  In 
order  to  estimate  the  effects  of  radiation  damage  the 
measurements  were  extended  to  investigate  the  temperature 
dependence  of  the  interaction. 

Two  experiments  in  addition  to  the  interaction  of 
Ga  in  Zn,  were  begun.  One  was  an  attempt  to  measure  the 
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quadrupole  interaction  of  Ga  in  Ga  metal,  the  other  an 
experiment  to  compare  the  quadrupole  interactions  of 
1 8  F ( 1 1 2 1  keV  level)  and  19F  (197  keV,  Q  =  0.11  b).  Both 
these  experiments  need  only  the  final  data-taking,  but  this 
was  prevented  by  a  major  shutdown  of  the  beam-line. 

For  the  experiments  various  requirements  were  present, 
especially  for  the  target  chamber.  The  Ga  in  Zn  experiment 
had  no  particular  requirement  for  the  room  temperature  run, 
but  for  the  study  of  temperature  dependence  a  hot  target 
chamber  was  developed,  which  was  capable  of  keeping  the  target 
at  a  well  defined  temperature  as  high  as  600  K,  did  not 
produce  background  radiation  and  showed  isotropic  absorption. 
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The  annihilation  radiation  of  the  decay  of  ^^i/2=  min.) 
produced  so  much  background,  that  it  was  necessary  to  shield 
a  part  of  this  radiation  from  the  detectors.  This  was  done 
with  the  background  supressing  target  holder,  which  uses  up 
to  5  targets  in  the  same  measurement.  The  Ga  in  Ga  metal 
experiment  required  a  cooling  device,  because  the  melting 
point  of  Ga  is  only  29.8°C  and  the  purpose  was  to  look  at  a 
solid  target.  It  also  required  the  use  of  a  less  efficient 
detector,  for  the  yield  of  the  reaction  of  interest  was  very 
high.  The  Low  Energy  Photon  Spectrometer  (LEPS)  detector 
was  used  for  this  purpose. 

All  experiments  were  analyzed  with  the  assumption  that 
the  interaction  originated  from  a  static  electric  field 
gradient.  This  assumption  implies  that  the  radiation  damage 
either  is  neglibible  or  is  "frozen  in"  for  a  time  long 
compared  to  the  lifetime  of  the  level  of  interest.  This  was 
justified  by  the  measurements  of  Bleck  et  al .  (B172)  and 
McDonald  et  al.  ( Me 7  2).  But  the  results  of  this  work 
indicate  that  the  field  gradient  experienced  by  the  nuclei 
is  probably  not  static  for  all  nuclei. 
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Fig.  1: 


Euler  angles  a  rotation  around  the  z  -  axis 
$  rotation  around  the  y-j-  axis 
y  rotation  around  the  z'-  axis 
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where  D^,  (a$Y)  is  the  rotation  matrix. 

When  we  look  at  systems  with  cylindrical  symmetry  around 
the  Z-axis  a  rotation  around  the  z-axis  will  not  change  the 
density  matrix.  Such  a  rotation  is  a  M»  M  0,  y  =  0  and 
DqXq,  (a, 0,0)  -  e-ic<a6qq, 

and 


(p')g,  =  E.  e'laq6 

H 


But 


Cp  1 ) 


q  -  qq '  rq 

-  i  q  '  a  A  x 
e  Pn.<5 

A 


qq' 


which  means  that  only  the  terms  with  q  =  0  do  not  disappear. 

The  remaining  terms  are  simply  related  to  the 
alignment  parameters  B  (I)  as  defined  by  Ro67  by 


(I) 


1 

(21+1 


Ba  (I) 


(II. 4) 


When  a  state  as  described  above  decays  with  a 
Y  -  ray  to  a  state  with  spin  we  have  an  angular  distribution 
as  derived  in  Ro67 


W(e  )  ~  £  ^  even  ^k  ^1^  R  k  ^1^2^  R  k  ^ ^ ® ^  ®  ^  ^ ) 

where  (cos  0)  is  the  Legendre  polynomial  in  the 
cosine  of  the  angle  with  the  z-axis. 

Rk  ( 1 1  1 2  ^  transition  dependent  parameter. 
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( 1 1  1  2  ^ 


=  E 


(  Ltt  )  (L  '  7T  '  ) 


x  <  TT  >  x  <  TT  '  > 

tRkUL,i1I2)  L  l'  1 


\\  «L<7r>  I 


(II. 6) 


Rk  (LL'I-jI^)  is  a  product  of  C 1  ebsc h-Gordan  and  Racah 
coef f i c i ents ( see  formula  3.31  of  Ro67). 

L  and  L‘  are  the  allowed  multipolarities  of  the  decay  and 
<tt>  indicate  an  Electric  or  Magnetic  transition. 


.  <  7T  > 


is  the  so-called  mixing-ratio  defined  by 


<7r>  _  <1,11  Tl<t.>  1  |  I2  >/(2L  +  1) 
L  <1,11  T-l<tt>  |  |I2  >/(2L  +  1) 


(II. 7) 


where  Tl<tt>  is  the  transition-operator  for  the  !  trans¬ 
ition  and  I  the  lowest-order  multipolarity  in  the  transition. 
The  angular  distribution  of  the  second  y  -ray  from  the  state 
with  spin  I2  to  the  state  with  spin  is  given  by 


W(0)  =  B  k  (.  I )  Rjl,!,)  Ub(I,I,)  P^lcos  6) 


k  '  2  3  '  k '  1  2 


(II. 8) 


wi 


th  Uk^Il12)  ZL,2  (6L,2)2  Uk^L12IlI2^  /  X<6L,2) 


Uk^L12IlI2^ 


k  W(I,I,I2I2;k  L,2) 
W(I,I,I2I2;0  L,2)  * 


(II. 9) 


where  W  is  the  Racah  coefficient  and 

l_-|2  is  the  multipolarity  of  the  y-decay  from  1^  to  1 2 
Usually  and  also  here  we  take  only  mixing  between  two  multi¬ 
polarities,  for  the  transition  probability  becomes  negligibly 
smaller  with  higher  multipolarity. 
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Formula  II. 8  is  only  valid  if  the  state  with  spin 
I2  is  unperturbed,  i.e.  the  lifetime  is  short  compared  to 
the  time  in  which  an  interaction  can  disturb  the  density 
matrix . 


For  practical  purposes  the  angular  distributions 
are  described  by 

W(9)  “  Ao  Zk  even  ak  Pk  <cos0)  <n-10> 


where 


VkRk 

BUR 

OOO 


(II. 11) 


(U^  =  1  when  we  observe  the  first  decay  and  is  the  value 
for  the  observed  decay). 

When  level  1  is  perturbed,  i.e.  has  a  time-dependent 
density  matrix,  we  end  up  with  a  time-dependent  angular 
distribution  for  both  decays. 


w(0)  =  A0  Ik  a k ( t )  Pk  (cos  0) 

=  £k  B k ( I ,  t )  Uk  R k  Pk(cos  6) 

=  Zk  B k( 1 , 0 )  Uk  R k  P k ( c 0 s  6)  Gk(t) 

=  Aq  £  k  a  k  G  k ( t )  U  k  R  k  P  k ( c  0  s  6)  (11.12) 

ak(t) 

Where  Gk(t)  5  i“(OT 
•  0 

From  this  we  see  that  the  second  decay,  for  which  in  general 
Uk  f  1,  has  the  same  time-dependence  as  the  first  decay,  only 


' 


the  unperturbed  value  a  k ( 0 )  will  be  different. 
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II . 2  Perturbed  Angular  Distribution 

The  expression  for  Gk(t)  can  be  evaluated  for  any 
case  with  a  known  interaction-Hamiltonian,  which  is  time- 
dependent  in  one  system.  From  (II. 2)  and  (II. 4)  follows 


Gk«t> 


ak^  =  BkU.t)  = 
a  k  ( 0 )  B  k ( 1 , 0 )  pk(0) 


(11.13) 


In  solving  the  expression  for  G  we  have  to  work  in  three 
coordinate  systems: 

(a)  The  laboratory  system,  which  is  defined  by 
the  plane  of  detection.  In  this  system,  ak(0)  and  G^(t)  are 
defined. 


(b)  The  "creation"  system.  In  this  system  we  know 
the  density  matrix  as  formed  by  the  beam  and  subsequent 
reaction,  or  any  other  method. 

(c)  The  crystal  system,  which  gives  us  the  inter¬ 
action-Hamiltonian,  defined  by  an  externally  applied  field 


or  hyperfine  fields. 

To  calculate  G.(t)  we  follow  the  following  procedure: 

K 

(1)  Transform  p*(t)  from  the  laboratory  system  to 
the  crystal  system. 

(2)  Form  the  density  matrix  in  the  crystal  system. 

(3)  Calculate  the  time-dependence  in  this  system  of 
the  density  matrix  as  function  of  the  density  matrix  at 


time  zero. 
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(4)  Expand  the  time  zero  density  matrix  in 
statistical  tensors. 

(5)  Transform  back  to  the  laboratory  system  and 
we  end  up  with 

pXq{t)  =  VX'  GiS'(t)  Vl0)- 

The  time-dependence  of  the  density  matrix  is  given  by  the 
equation  of  motion. 

ihp(t)  =  [H,  p] 

We  consider  only  time- i ndependent  Hamiltonians  so 

p(t)  =  e'iHt/fi  p(0)  eiHt/R  (11.14) 

We  also  use  the  formulas  11.1,2,3. 

Pq  !  ( t )  =  Zq,  ,  (p  1  )q'.  ,  (t)  D^!  ,ql  (U) 

=  Eq,,m(I=m  I  m '  |X'q")(-l)A''1'm 

<1  m  |  p  •  ( t )  1  I  m  1  >  !  i  ,  (£2 ) 

=  £  i  .  ii  ii,  ( I-m  I  m '  | A  '  q  '  '  ) 
q ' 'm  m  '  '  1  ^  ' 

(-l)X''I'm  D  X  .  q ,  ( S2 )  x 

x<I  m  |  e*iHt/h|I  m  1  1  ><  I  m  '  '  |  P  '  ( 0 )  |  I  m1'^  x 

x  <1  m"  '  |eiHt/R|l  ■•> 


. 


•:)'->  .'\a 


( 
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Zq"  mm'  'm  '  '  1  A  ^ 1  "m  1  m '  I  A'  9  "  )  * 

X  (-l)A_1"m  DA',q,  (n)x 

x  <1  m  |  e~lHt/^  |I  m"  ><  im 1 1 '  |elHt/^  |I  m‘>  x 
x  (p')qll  ,(0)  (I  -m"  I  m"  '  |Xq"  ,)(-l)X'I'm" 


v  fTmTm*  i  a  i  ^,11  \  /  t  sA+A-2 1  -m-m 

V'qmm"m"  '  X(I'm  1  m  lx  P  )(-l)  x 

x  <1  |I  m">  <1  m"  '|  eiHt/fi  |I  m'>  x 


x  (I  -m"  I  m"  '|X  q"  •)  pA(0)  DqAql ,  .  (-«)  DA ql  (£2) 


h  IAq  gJJ,'  (t,Q)  pj(0) 


(11.15) 


This  is  the  general  expression  for  the  pertubation  coefficient 
due  to  a  static  Hamiltonian.  It  is  applicable  in  all  cases 
of  geometry. 

There  are  a  couple  of  experimental  conditions,  which 
can  reduce  this  expression  to  more  workable  formulae. 

(a)  The  'creation'  system,  in  which  p  is  cyl i ndr i ca 1 1 y 
symmetric,  coincides  with  the  laboratory  system.  That  means 
pq  ( 0  )  =  0  for  q  f  0 . 


So 


Pq.(t)  =  ^ 


x  %(0> 


(11.16) 


i 


1  3 


(b)  A  random  orientation  of  the  crystal  systems  with 
respect  to  the  laboratory  system,  as  encountered  in  micro- 
cystalline  samples.  Then  we  can  perform  an  integration  over  Q 


1 

4  7T 


/  D 


oq 


I  I 


(-n)  d  , .  ,  (n)  dsr2  = 


i 

4  7T 


*  \ 

/  D  , , 

q 


.0(n)  d  m  ,  (a)  dft  = 


l 

2A  +  1 


<5  , 

q  o 


(H.17) 


Here  we  can  integrate  over  two  angles  a  and  3,  for  y  may  be 
taken  zero  and  still  we  will  describe  the  rotation.  When  we 
are  averaging  over  the  total  solid  angle  the  factor  4 tt  drops 
out.  Substitute  (11.17)  and  (11.16)  in  (11.15) 


p*u) 


„  /T  T  i  | ,  1 1  \  /  *i  X2A  -  21  -  m  -  m'  ' 

=  Eq'-mn"  (J  "m  1  m  lX  <1  )(“!)  x 

x  <1  m|  e"lHt/^  |I  m">  <1  m"  '  |  elHt/^  |I  m’>  x 

(I  -m"  I  m"  '  | X  q"  )  p*(0)  2X\  ^  (11.18) 


(c)  The  interaction  Hamiltonian  is  diagonal  in 


|  I  m> 


,  I  —  1  H  t/ fi  I,  m,  ,  .  .  _  i  H  t/  R  t 

<  I  m  |  e  |  I  m  >  =  e  6 


mm 
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This  substituting  in  (11.18) 


po(t)  =  VJ1  'm  1  m'|x 


2A-2I  -2m  -iE  t/h 
e  m 


x  elEm,t/h  (I  -m  I  m‘|Xq")  pA(0)  1 


Sq.  .m  (I  -m  I  m  •  |  q")2cos{(Em  -  Em ,  )  t/R}  pA(0) 


(11.19) 


^  1  j  2  X  -  2 1  —  2m 


=  1  for  the  exponent  is  always  even. 


(Em  ~  Em ' ^  =  ~(Em'  "  Em ^  so  the  sine  functions  will 
cancel  each  other  and  we  are  left  with  the  cosines. 

Formula  (11.19)  is  the  often  derived  formula  useful 
for  an  axially  symmetric  quadrupole  interaction,  measured  in 
a  microcrystalline  sample.  For  a  magnetic  field  (Em  -  Em , ) 
is  a  constant  for  the  allowed  transitions  of  Am  =  ±1,  and 
this  will  further  reduce  the  expression. 

For  a  Hamiltonian,  which  is  not  diagonal  in  the  |I  m> 
representation,  we  have  to  diagonalize  it  to  a  set  of  eigen¬ 
vectors  |n>.  The  transformation  from  |I  m>  to  | n>  is  described 
by  the  transformation  matrix  U.  U  is  unitary  and  U  =  1. 


rrm 


<m|  H  |m'>= 


Z  ,  <m|  U_1U  H  U_1U  |m '>= 
mm'  1  1 


^  m' ' m1 


nn  mm 


<m'  '  |  n>  <m ' '  |  m ' ' >  <n|  H  |  n ' ><m' '  '  |  m"  ' ><n '  |m‘>  = 


* 
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Z  i  U  E  6  ,  U  , 

mm  n  mn  n  nn  n  m 


Z  ,  U 
mm  n  mn 


E  6  ,  U*. 

n  nn'  m  n 


Use  in  (11.18) 


q*  1  mm' '  nn' 


(I  -m  I  m‘  |A  q")  (-1)2X  21  *  m  ‘m 


x(I  -m"  I  m"  '  |X  q"  )  p*(0) 


(11.20) 


This  is  the  formula  for  the  perturbed  angular  distribution 
for  a  non-axially  symmetric  Hamiltonian,  from  a  cyl i ndr i ca 1 1 y 
symmetric  state  at  time  zero  i n a  microcrystal  1 i ne  sample. 

As  is  evident  the  pertubation  factor  G(t)  =  p(t)/  p(0)  is 
only  dependent  on  the  multipolarity  of  the  y-decay  and  the 
Hami 1  ton  i  an . 

1 1 . 3  The  Hamiltonian  for  Quadrupole 
Interaction 


The  electrostatic  Hamiltonian  of  a  charge  distri¬ 


bution  in  an  external  field  is  (De53) 

Hel  =  /  p(x)  V  (x)  d3x 

Where  p(x)  the  local  charge  density  is. 

A  Taylor  expansion  of  the  potential  about  the  origin  is 


V  (x )  =  V(0)  +  x.VV(O)  +  Ij 


(0)  xiXj  + 


U  .  "  Z 


' 


Vv  I 


So  Hgl  =  /  p(x)  V  (0)  d3x  + 

+  /  p(x)  x.  ^V(O)  d3x  + 
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+  /  P(x)  xiXj  (0)  d3x  +  . 

(11.21) 

The  first  term  is  the  energy  of  the  nucleus  in  its  surround¬ 
ings  and  not  important  for  this  problem.  The  second  term 
vanishes  because  the  nucleus  is  in  a  minimum  of  the  potential 
so  V  V  =  0.  The  third  tern  is  the  quadrupole  interaction,  the 
higher  order  terms  are  neglected.  We  can  write  the  quadrupole 
interaction  as 


H 


t!2> 
i  j  i  j 


(11.22) 


where  T^^  is  the  second-rank  tensor,  the  quadrupole  moment 
tensor,  and  T .  ^  =  f  p(x)  x^x^  d3x  and  =  ^x~^x~  (0)  is  the 

component  of  the  field  gradient  tensor. 

It  is  always  possible  to  make  the  field  gradient  tensor 
diagonal  by  a  suitable  choice  of  the  coordinate  system. 
Furthermore  the  Laplace  equation  is  valid  for  V  describes 
only  the  external  field. 


V  +  V  +  V  =  0 
xx  yy  zz 


(11.23) 


By  defining  the  axis  properly  the  usual  condition  is 


Define  n  =  (Vxx  -  V  )/V 

xx  yy"  zz 

From  (11.23)  and  (11.24)  it  follows  that  0  <n<  1. 

way  the  field  gradient  is  described  by  Vzz  and  n  and 

orientation  of  the  coordinate  system. 

(2)  (2) 

Tv  '  and  Vv  '  can  be  expanded  in  spherical  tensors 

t<2>  =  I  +  2 

q  =  -2 

\/(2)  =  i  v2 
q  q 

SO  H  -  i  Eq  T2 

For  a  diagonal  field  tensor 

v±i  =  0 

V0  -  Vzz 

V2  =  —7—  (V  “V  )  =  )  r  n  V 

v±2  /6  v  XX  yy'  76  zz 

<1  m |  H  | I  m’>  = 

<1  m|i  Tq  V?q  I1  "  ’  >  = 

I  zq  V-q  <!  ml  Tq  I1  m ‘ >  = 

i  I  (-I)4  <I|  |  T(Z)  |  |I>  V2  (Im'  2  q|I  m 

q  • 

The  quadrupole  moment  of  the  nucleus  is  defined  as 
eQ  =  /  Pn  (r2  - 


and 

q 
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(11.25) 

In  this 
the 


(11.26) 


z2 )  d3x 


■' 


18 


=  2  /  Pn  x2  d3x 

=  2  <11 |  Tq  j  1 1 > 

=  2  (-1  )"  (I  I  2  0 1  I  I )  <  I  |  |  T^2M|I> 

So  <Im|H|lm'>= 

i  2q  eQ  (I  m '  2  q | I  m)  (I  I  2  0 | I  I)-1  V^q 

Calculate  the  Cl ebsch-Gordan  coefficients  and  substitute  the 
values  for  the  components  of  V 

q  =  °  <  I  m  I  H  1 1  m '  >  »  (3m2  -  1(1  +  1))  _V)-  Vzz  6mm, 

q  =  ±1  <1  m|  H  | I  m‘>  =  0  (11.27) 

q  =  ±2  <Im|H|Im+2>= 

{(I  +  m  -1)(I+  m )  ( I  ±  m  +1)(I±  m+2)}*  ^  V  zz 

These  values  are  also  derived  by  Matthias,  Schneider  and 
Steffen  (Ma62  and  Ma63)  using  the  expansion  of  the  Coulomb 
potential  in  spherical  harmonics.  It  should  be  noted  that 
the  definitions  in  this  reference  for  both  T  and  V  differ  a 
constant  factor  from  the  definition  used  above. 
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CHAPTER  III 


EQUIPMENT 

Beam 

We  used  the  7  MeV  Van  de  Graaff  accelerator  of  the 
University  of  Alberta,  with  a  pulsed  beam  of  1  ns  width 
with  a  repetition  rate  of  1  MHz. 

De tec  tors 

As  detectors  we  used  4  Ge(Li)  detectors,  2  in  each 
experiment,  with  volumes  of  20.7,  23.0  and  27.4  cc  and  the 
LEPS  detector.  The  three  large  detectors  combine  good 
energy  resolution  with  good  time  resolution  and  a  wide 
efficiency  range.  The  LEPS  detector  (3  cc)  is  very  good  for 
low  energies  (less  than  400  KeV )  and  less  sensitive  to  higher 
energy  y-rays  reducing  the  c ou nt i ng -ra te  by  a  large  factor. 

Electronics 

We  had  two  electronic  setups,  one  for  two  time 
dependent  measurements  at  one  time  and  one  for  one  timing 
and  one  monitor.  The  double  timing  had  the  drawback  that 
the  tuning  was  more  complicated  and  that  the  sum  boxes  were 
not  able  to  handle  high  counting  rates.  The  single  timing 
cost  more  time,  but  was  easy  to  set  up,  and  could  handle  high 
counting  rates.  The  time-resolution  of  the  system  was 
typically  2.7  ns  FWHM  for  843  keV  and  8.0  ns  for  188  keV . 
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The  electronic  arrangement  with  two  timing  counters 
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Fig  3:  The  electronic  arrangement  with  the  single  timing  counter 
and  the  monitor  counter 
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A  major  problem  which  showed  up  was  the  appearance 
of  10  and  30  MHz  ripples  in  phase  with  the  1  MHz  beam  pick-off. 
This  showed  up  in  the  time  spectra  as  a  modulation  of  the 
yield  as  function  of  the  time.  We  suppressed  it  partly  with 
a  50  ns  delayed  reflection,  gaining  at  least  a  factor  4. 

For  the  timer  we  used  a  time- to-ampl i tude  converter 
(TAC),  supplied  with  a  start-signal  from  the  Ge(Li)  detector. 
The  stopping  side  was  from  the  beam  pick-off.  This  signal 
was  delayed  to  place  the  prompt  peak  on  the  end  of  the 
spectrum.  This  set-up  had  the  effect  that  there  was  always  a 
stop-signal  and  that  the  time  in  our  time  spectrum  was 
running  backwards. 

The  calibration  was  performed  with  a  commercial 
calibrator  (Ortec  model  462).  For  on-line  sorting  of  data 
we  used  the  program  TWO  PARAMETER  SORTING,  available  on  the 
Honeywell  516  computer.  For  description  see  Chapter  IV. 

Target  Mounting 

(1)  Normal  beamstop  target.  In  the  figure  4  are 
also  indicated  the  other  devices,  i.e.  pulse  pick-off  which 
gives  the  time-signal  from  the  beam,  the  energy  slits  and  the 
collimator,  together  with  the  definition  of  0  .  The  other 
target  holder  assemblies  replace  that  shown  to  the  right  of 
the  dotted  line  in  Figure  4.  The  collimator  is  not  shown 

in  the  next  figures 

(2)  Hot  target  chamber.  For  heating  up  experiments 
we  designed  a  target  chamber  with  a  heating  system  and  control 


' 
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Fig  4:  The  normal  beam  stop  target  holder  with  surrounding 
devices 
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unit.  After  some  trials  this  final  set-up  had  these  properties: 

(a)  Ceramics  were  high  temperature  castable  ceramics 
from  Aramco. 

(b)  The  target  holder  can  be  placed  at  an  angle  up 

to  30°  to  the  beam,  made  possible  with  a  rotatable 
connection  above  the  hea ti ng -wi r e . 

(c)  The  target  holder,  the  heating  wire  and  the 
nickel  hat,  which  served  as  a  Faraday  cup  and  as 
a  shielding  of  the  surroundings,  were  electri¬ 
cally  isolated  to  the  ground  and  each  other. 

(d)  The  obtainable  temperature  was  at  least  600  K. 

(e)  The  required  power  input  for  560  K  was  approximately 
35  Watt. 

(f)  The  temperature  was  stable  within  0.1  K,  for  each 
temperature,  if  the  beam  was  constant.  Large 
fluctuations  of  the  beam  resulted  in  the 
temporary  change  in  temperature.  Too  much  power 
input  resulted  inlarge  temperature  variations. 

(g)  Due  to  the  heat  input  from  the  beam  the  temperature 
of  the  centre  of  the  target  was  70  -  100  K  above 
the  measured  temperature.  This  is  calculated 

in  Appendix  I  and  verified  with  the  sublimation 
rate  of  natural  zinc.  All  quoted  temperatures 
are  85  K  above  the  measured  temperatures. 

(h)  The  thermocouple  was  a  standard  Fe/Co ns ta nta n 
thermocouple.  The  results  were  interpreted  with 
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Fig  5:  Hot  target  chamber.  The  5  feedthroughs  are  for  the  heating-wire 
(2)  the  thermocouple  (2)  and  the  current  collector 
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Fig  6:  The  control  box  circuit  diagram.  This  circuit  switches  the 
power  on  or  off  if  the  reading  of  the  thermocouple  is  lower 
or  higher  than  required 
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the  tables  in  the  Handbook  of  Chemistry  and 
Physics  (Ha72). 

(3)  Background  surpressing  target  holder.  For  the 

1  8 

F  experiments  we  encountered  a  high  511  keV  y-yield  from 

+  1 8 

the  annihilation  of  the  3  of  the  F  decay.  Using  the  fact 

that  the  half-life  is  110  min,  we  shielded  in  this  set-up  four 

targets  with  2  cm  lead  and  used  the  fifth  target  for  the 

measurements.  By  changing  the  targets  every  30  min.  we  were 

able  to  reduce  the  compton  background  from  the  511  keV  y  ray 

in  the  region  of  interest,  180  -  200  keV,  by  a  factor  3  to  4 

1  8 

relative  to  the  yield  of  the  188  keV  F  y  ray.  The  ladder 
with  the  five  targets  can  slide  in  and  out  without  disrupting 
the  vacuum.  The  targets  are  electrically  isolated  from  the 


chamber . 


.. 
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Fig.  8:  The  background  suppressing  target  holder. 
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CHAPTER  IV 


DATA  ANALYSIS 

As  described  in  Chapter  III  the  electronics  give  to 
the  Honeywell  516  computer  3  signals,  an  analog  time-signal, 
an  analog  energy-signal  and  a  logic  'event' -pulse.  The 
event-pulse  serves  the  purpose,  that  during  the  pulse  the 
computer  accepts  both  other  pulses  and  processes  them 
according  to  the  program.  ADC's  convert  both  other  signals 
from  analog  to  digital  with  up  to  4096  channels  after  which 
the  program  TWOP  stores  both  in  a  separate  memory  region  and 
if  the  energy  happens  to  fall  in  one  of  the  set  regions  in 
the  energy  spectrum,  the  so-called  windows,  its  corresponding 
time  is  also  stored  in  a  third  memory  region.  So  we  have 
three  spectra,  the  energy,  the  times  of  all  y  ray's  and  the 
times  of y  ray's  with  certain  energies,  for  each  window  one 
complete  time-spectrum.  After  each  run  the  data  was  written 
on  magnetic  tape  for  off-line  analysis  on  the  XDS  940. 

The  off-line  analysis  is  with  the  program  TDPAD 
(Lu75)  especially  developed  for  the  analysis  of  these 
experiments.  It  is  able  to  fit  time  spectra  and 
background  of  the  peak  of  interest  to  the  time-dependent 
angular  distribution  coefficients  and  fits  that  to  the 
theoretical  curve  which  depends  on  the  f requency , . i ts 
relative  spread,  the  spin,  the  asymmetry  parameter  n»  the 
unperturbed  value  of  a^  and  a^,  the  AQ(t)  and  the  interaction 
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frequency  w  =  (e2Qq)/h. 

The  TDPAD-program  is  essentially  fitting  to  the  formulae 
(II. 11)  and  (11.20),  using  the  Hamiltonian  given  by  (11.27). 

The  pertubation-coeff icient  as  described  is 

Gk(t)  =  EN  skN  C0S  (%*)  (IV. 1  ) 

For  the  axia I ly- symmetric  case  N  is  integer,  but  deviates 
in  general  from  integral  values  for  n  t  0.  The  s.  are  easily 
obtainable  from  (11.20), 

N  =  (E  -  E  ,)/(ho)  )  and  oj  =  (e2Qq)/h. 
n  n  q  q 

Because  all  nuclei  do  not  feel  the  same  field  gradient, 
due  to  slight  displacements  in  the  microcrystals,  we  have  a 
relative  frequency  spread  6w.  We  assume  the  ^Lorentzian 
frequency  distribution  given  in  (Ha73). 

Gk(t)  =  EN  skN  e"N(°q  6l0t  cos  NlV  (IV-2) 

Both  (IV. 1)  and  (IV. 2)  are  the  theoretical  curves  for  infinite 
time  resolution.  This  is  experimentally  unobtainable,  so  we 
have  to  take  into  account  the  time  resolution  of  the  system. 

This  is  a  Gaussian  distribution  and  can  be  described  by 

Gk'(t)  =  EN  skN  e'N“q6wt  e(NV)2/2  (cos  Nuqt) 

(IV. 3) 


where  t  is  the  time  resolution. 
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The  fit  is  to 


ak(t)  =  ak(0)  Gk'(t)  +  c,  where  c  is  a  time- 

independent  constant.  Gk(t)  is  calculated  by  numerically 
diagonalizing  of  the  Hamiltonian. 

The  theoretical  curves,  as  they  are  calcualted  with 
the  formulas  (IV. 1),  (IV. 2)  and  (IV. 3)  are  shown  in  the 
figures,  for  a  couple  of  typical  spins  and  spins  of  interest, 
with  various  asymmetry  parameters,  resolution  times  and 
frequency  spreads. 

In  the  literature  is  often  given  instead  of 
is  noted,  w  is  defined  as  the  smallest  energy-difference 
for  an  axially  symmetric  field 


“o  =  41(21- 1  )'  wq 
“o  =  41(21-1  )  “q 


for  integer  spin 


for  half-integer  spin. 


We  prefer  to  use  w 
action  parameters 


q 

Q 


for  it  is  dependent  only  on  the  inter- 
and  q  and  not  on  the  spin. 


-3  Iti  6  i  9^  5  ♦  (  ■  '  (0  \}6  s  iB  ' 
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Figures  9  -  18: 

The  values  of  the  G2(t)  term  as  function  of  time  for 
different  values  of  spin,  n»  Sw  and  resolution  time 
t.  Note  that  Fig.  9  is  the  G ^ ( t )  curve  for  spin  2. 
For  all  curves  the  first  case  is  represented  by  the 
solid  curve,  the  second  by  the  dashed  and  the  third 
by  the  dash  dot  line. 
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Fig.  12:  Spin  =  5/2  n  =  0.0, 0.5, 1.0  6w  -  0.0  x  -  0.0 
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Fig  15:  Spin  =  4  n  =  0.4  6<d  =  0.0  t  =  0.0, 1.0, 2. 5 

Channels  for  a)  =  5  rad/ch. 


Fig  16:  Spin  =  4  n=0.0  6u)  -  0.0  t  -  0.0, 1.0, 2. 5 

Channels  for  u>  =  5  rad/ch. 
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Fig  17:  Spin  =  4  n  =  0.4 


6w  =  0,2,5%  x  =  0.0 


Fig  18:  Spin  =  4 


n  =  0.0 


6oj  =  0,2,5% 


x  =  0.0 


CHAPTER  V 


THE  EXPERIMENTS 

7  0  7ft 

V  .  1  Zn  ( p  ,  n )  Ga 

The  experiment  was  initiated  because  a  discrepancy  existed 
between  the  spin  assignment  from  particle  work  and  the  measured 
anisotropy  in  the  y  decay  of  the  879  keV  level,  especially 
when  compared  with  the  411  keV y  ray  from  the  1102  keV  level, 
which  also  was  assigned  a  spin  4  (Do72).  Together  with  the 
life  time  and  the  possibility  of  an  Electric  Field 
Gradient,  we  extracted  the  unperturbed  value  of  the 
anisotropy,  to  see  if  this  is  in  agreement  with  the 
prediction. 

The  879  keV  level  of  Ga  has  a  half-life  of  22.7 
ns  ( H  u  7  5 )  and  decays  to  the  691  keV  level  (I71  =  2  ).  Both 
the  188  keV  and  the  691  keV  y  rays  are  suitable  for  doing  a 
time-dependent  measurement  as  long  as  one  takes  into  account 
that  the  region  of  188  keV  is  populated  with  much  abundant 

19  197  64 

y's,  i.e.  197  keV  of  F,  191  keV  of  Au ,  185  keV  of  Cu, 

7  1 

175  keV  of  Ge  and  the  bac ksca tter i ng  edge  of  the  Comptons 
of  511  keV.  The  691  keV  coincides  with  the  692  keV  E0 
transition  in  Ge  from  (n,n')  in  the  Ge(Li)  detectors. 

Another  drawback  is  the  relatively  short  lifetime  of  the  decay. 
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Fig.  19:  The  level  scheme  of  70Ga  with  the  transitions  of 
i nterest . 
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Figure  20: 

A  typical  y  -  spectrum  from 


7  o 

(p,n)  Ga. 


labels  are  the  energies  in  keV. 


The 


counts 
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In  developing  a  hot-target  chamber  the  aluminum  has 
to  be  shielded  for  this  gives  too  much  background  radiation. 

The  target  should  also  be  prevented  from  evaporating.  The  use  of 
mica  for  this  purpose  is  unacceptable  for  it  gives  too  much 
annihilation  radiation.  The  solution  was  found  in  clamping 
two  .0003"  thick  nickel  foils  around  the  target. 

7  0 

The  target  was  a  1.9  mg/cm2  foil  of  enriched  Zn 
(68%).  The  beam  was  2.85  or  2.9  MeV  protons,  which  is  300 
keV  above  the  threshold  for  the  879  keV  level,  and  the  nuclei 
are  aligned  by  the  reaction.  There  were  three  measurements 
done  in  a  single  run,  one  at  410  K,  one  at  420  K  and  one  at 
510  K.  We  used  the  doubl e- timing  system  and  measured  only  two 
angles,  0°  and  90°.  This  has  the  problem  that  the  &2  Term  is 
determined  only  when  a^  term  is  negligible. 

A  run  at  560  K  and  one  at  340  K,  were  performed 
using  ..the  single  timing  set-up  and  at  five  angles  each. 

We  fitted  the  data  with  a  normalization  to  the  651 
keV  line,  assuming  that  the  ti me- i n tegra ted  yield  of  this  decay 
is  isotropic  ( Do  7  2 ) .  The  time  zero  is  defined  by  the 
position  of  the  background  prompt  peak  in  the  time-spectrum, 
fitted  with  peakfit  program.  The  time  resolution  was  measured 
with  a  prompt  peak  in  the  region  of  interest. 

When  the  data  were  analyzed  with  an  infinite  time- 
resolution,  they  yield  results  published  in  Nucl.  Physics 
(Hu75),  and  reproduced  in  Table  1.  These  values. have  one 
feature  which  makes  them  unbelievable,  the  high  frequency 
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UNPERTURBED  VALUES  OF  LEGENDRE  POLYNOMIAL  COEFFICIENTS 
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spread.  When  we  introduce  the  finite  time  resolution  and 
the  asymmetric  field  gradient,  it  turns  out  that  the  time 
resolution  pushes  the  frequency  spread  down  to  zero. 

Out  of  the  analysis  it  followed  that  one  cannot  distin¬ 
guish  clearly  between  the  various  spins.  Neither  spin  3  nor 
spin  5  can  be  excluded  for  the  879  keV  level  with  these 
experiments.  This  is  caused  by  the  short  life  time  compared 
to  the  full  period  of  the  time-dependence.  Only  the  first 
recovery  of  the  anisotropy  is  seen  and  this  is  not  sufficient 
to  determine  the  spin.  Only  spin  3,  4,  5  were  consider ed  for 
the  others  are  excluded  on  the  basis  of  spin  assignments 
from  particle  reactions.  All  the  values  calculated  below 
are  for  spin  4.  If  spin  3  or  5  is  chosen  for  the  level,  one 
must  scale  these  values.  When  the  recovery  is  set  at  a  fixed 
time,  o>  will  be  different  for  the  three  cases.  Then  o)Q  = 

3w  / (41 ( 21 -1 ) ) ,  so  a)  is  certainly  different  for  the  three  cases. 
•  q 

TABLE  2 

SCALING  FACTORS  FOR  u  FOR  DIFFERENT  SPINS 


Spin  a)  oj 

r  o  q 


3  1.520  0.814 

4  1.000  1.000 

5  0.847  1 .361 

(n  =  0,  values  normalized  to  spin  4) 


■ 
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This  is  confirmed  with  the  analysis  of  the  data  with  the 
three  different  spin  values. 

Neither  is  the  value  of  n  determined.  For  spin  4  the 
value  of  n  is  restricted  to  0.0  <  n  <  0.6.  This  large  range 
is  also  due  to  the  small  part  of  the  full  cycle  which  is 
visible. 

All  values  of  a^  and  a^  can  be  affected  by  a  mis¬ 
alignment  of  the  set-up  up  to  0.05.  This  error  is  not  quoted 
in  the  table,  but  is  common  to  the  188  and  the  691,  and  will 
appear  in  the  constant  term  c,  which  is  a  time- i ndependent 
constant  added  to  a^G^t)  to  get  "a£(t)".  This  constant  term 
absorbs  also  the  normalization  errors. 

All  data  or  typical  samples  of  it  are  shown  in  the 
table  and  figures.  The  a^  of  the  691  y  ray  is  consistent 
with  0.00  and  not  quoted  in  the  table. 

V . 2  The  F  Trial 

We  tried  to  measure  the  quadrupole  moment  of  the  1121 

keV  level  of  ^F,  which  decays  with  a  184  keV  y  ray  to  the 

937  keV  level,  which  subsequently  decays  to  the  ground  state. 

The  1121  keV  level  has  a  spin  5  and  a  half-life  of  145  ns. 

Since  fluorine  doesn't  have  a  simple  non-cubic  crystal 

1  8 

in  any  of  its  compounds,  we  wanted  to  compare  F  with  the  197 

1  9 

keV  level  in  F,  which  has  a  known  Quadrupole  moment. 
Radiation  damage,  which  is  dependent  on  the  beam  intensity  and 
energy,  can  change  the  effective  E.F.G..  So  we  wanted  to 


, 
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Table  3 


THE 

DATA  FROM 

THE 

VARIOUS 

RUNS  FOR  7  °Ga 

188  keV 

ray 

T(K) 

o)q(MHz) 

X2 

Nf 

C 

a2(0) 

a2t0) 

560 

330 

63.7 

34 

0.04 

0.278  ±  0.008 

0.32  ±  0.01 

510 

316 

41.5 

25 

0.13 

0.255  ±  0.026 

0.38  ±  0.03 

420 

315 

28.2 

25 

0.17 

0.264  ±  0.021 

0.43  ±  0.02 

410 

316 

28.7 

25 

0.17 

0.209  ±  0.020 

0.38  ±  0.02 

340 

319 

19.5 

24 

0.08 

0.341  ±  0.022 

0.42  ±  0.02 

691  KeV 

ray 

560 

330 

25.6 

23 

-0.13 

-0.156  ±  0.017 

-0.29  ±  0.02 

510 

316 

37.4 

26 

-0.06 

-0.169  ±  0.020 

-0.23  ±  0.02 

420 

315 

32.9 

26 

-0.05 

-0.170  ±  0.018 

-0.22  ±  0.02 

410 

316 

17.6 

26 

-0.03 

-0.186  ±  0.020 

-0.22  ±  0.02 

340 

319 

14.7 

22 

-0.025 

-0.257  ±  0.022 

-0.29  ±  0.02 

188  KeV 

ray 

a4(°) 

ai(0) 

560 

330 

41.2 

34 

-0.04 

-0.051  ±  0.014 

-0.09  ±  0.02 

340 

319 

12.1 

24 

0.00 

-0.144  ±  0.050 

-0.14  ±  0.05 

Nf  is  the  number  of  degrees  of  freedom  in  fitting  a(t)  to  the  theory 
a  1 ( 0 )  =  a(0)  +  C 

Time-integrated  values  for  the  188  keV  transition. 

Temperature  ^  a4 

0.20  ±  0.01 
0.21  ±  0.01 


340  K 
560  K 


-0.11  ±  0.02 

-0.06  ±  0.02 


, 

Fig.  21-34: 


The  experimental  results  of  ^UGa  in  Zinc  at 
various  temperatures. 
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a2G2(t) 


340  K  188  keV 
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340  K  691  keV 


a2G2(f ) 


t  (ns) 


a2^2(*' ) 


t(ns) 


410  K  188  keV 


53 


a2G2(t) 


a2G2(f) 


•■(ns) 


420  K  188  keV 
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510  K  188 keV 
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a4G4(t) 


■ 
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560 K  691 keV 
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measure  both  frequencies  in  the  same  run,  to  ensure  that 
both  E.F.G.  are  the  same. 

16  18  1  9 

We  used  the  reactions  0  (3He,p)  F  and  F  (3He,3He) 

1  9 

F*  with  2.95  MeV  3He  beam.  The  target  was  an  oxide  with  a 
layer  CaF^  evaporated  on  it.  The  present  recoil  energy 
required  a  thickness  of  CaF2  less  than  100  yg/cm2.  We 
measured  the  thickness  of  a  sample  with  a  particle  back- 
scattering  and  used  the  yield  of  the  197  keV  y  ray  of  this 
sample  as  calibration.  In  this  way  one  is  able  to  tell  if 
the  targets  are  thin  enough.  5  targets  are  less  than  100 
yg/cm2  in  thickness.  The  choice  of  the  oxide  is  bound  by 
some  conditions.  The  anion  in  the  oxide  should  not  give  a 
lot  of  radiation  in  the  3  MeV  3He  beam.  This  excludes  Si 
and  lower  Z  materials,  for  most  of  them  yield  8+-unstable 
nuclei,  producing  too  many  annihilation  y  rays  and  their 
comptons.  Secondly  the  oxide  should  be  non-cubic,  preferably 
hexagonal,  rhombic  or  tetragonal.  We  found  the  solution  in 
V205  molten  on  a  tantalum  backing,  with  CaF^  evaporated  on 
top . 

1  8 

The  F  annihilation  peak  produced  that  much  back¬ 
ground  in  the  180  -  200  keV  region,  that  a  background  surpass¬ 
ing  target  chamber  had  to  be  designed  using  the  fact  that  the 

1  8 

half-life  of  F  is  110  min.  One  target  was  used  for  30  min. 
then  shielded  for  the  detectors  while  measuring  on  the  other 
targets.  This  way  one  is  able  to  reduce  the  background  at 
184  keV  by  a  factor  3  to  4. 
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6  9  7  0 

V  .  3  Ga  ( d  ,  p ) _ Ga 

With  this  experiment  we  start  to  get  some  comparison 
between  gallium  ions  in  Zn  and  Ga  metallic  surroundings. 

This  already  has  been  done  with  Ge  as  probe  (Ha73a).  We 
used  natural  Ga  and  excited  the  879  keV  level  with  5.5  MeV 
deuterons.  The  only  experimental  difficulty  was  the  low 
melting  point  of  gallium  metal  which  is  only  29.8  °C.  The 
target  was  cooled  by  using  the  background  surpressing  target 
chamber  and  immersing  the  rod  of  the  target  ladder  into  ice 
water.  Probably  due  to  poor  thermal  contact  between  the 
target  and  the  ladder  we  didn't  succeed  in  maintaining  our 
target  solid,  and  measured  so  the  time  differential  pattern 

7  0 

of  Ga  in  liquid  Gallium.  The  results  show  the  feasibility 

7  0 

of  studying  the  quadrupole  interaction  of  Ga  with  the 
(d,p)  reaction.  The  unperturbed  value  of  a^  is  approximately 


0.1  . 
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CHAPTER  VI 


DISCUSSION  OF  THE  RESULTS 
V 1 . 1  The  Frequencies 

The  frequencies  quoted  in  Table  3  have  an  uncertainty 
of  3%,  in  which  region  the  X  changes  by  only  up  to  3%  of 
the  quoted  minimal  value.  A  change  in  the  frequency  will 
slightly  affect  the  unperturbed  a^  values. 

The  values  with  their  uncertainties  are  plotted  in 
Fig.  35  as  function  of  the  temperature.  From  this  graph 
follows  the  value  of  w  at  300  K  of  (315  ±  10)  MHz.  The 

q 

change  in  frequency  from  300  K  to  600  K  is  (+3  ±  6)%. 

A  definite  calculation  of  the  E.F.G.  has  not  yet 
been  performed.  According  to  the  current  theory 

eq  =  (1  -  >Jethat  +  (1  '  R)eqloc  (VI. D 

where  qlQC  is  due  to  the  electrons  in  a  shell  around  the 
nucleus  not  containing  other  ions, 

a  ^  is  due  to  the  other  ions  in  the  lattice  and  the 
Mlat 

electrons  outside  the  shell. 

represents  the  effect  of  the  displacement  of  the 
core  electrons  in  the  lattice  electric  field.  It  is 
called  the  Sternheimer  a nti -shi el di ng  factor. 


59 


' 


60 


T(K) 


Fig  35: 


The  dependence  of  w 


q 


on  the  temperature  of 


7  0 

Ga  in  Zn 
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R  is  the  shielding  factor,  which  expresses  the  change 
in  the  E.F.G.  of  the  electrons  at  the  nuclear  site 
due  to  the  interaction  of  the  nuclear  quadrupole 
moment  on  the  core  electrons. (Wa65  and  refs,  therein). 

Plat  can  be  calculated  for  various  lattices  ( W e 6 1 
and  Da6l),  yielding  for  our  hexagonal  Zn  lattice  eq,  .  = 

I  d  t 

1  6  2 

-1.5  x  10  V/ cm  .  y^,  the  anti -shi eld ing  factor  has  been 

7  0 

also  calculated  and  is  -9.905  for  Ga .  So  eq  -j  a  t  (1  -  ym)  = 

1 7  2 

-1.63  x  10  V/ cm  .  The  factor  R  is  calculated  by  Sternheimer 

for  various  ions  and  determined  to  be  0.1  to  0.4  (St63  and 

refs,  therein),  but  calculations  of  q^oc  have  not  yet  been 

performed  adequately  for  most  ions. 

For  an  estimate  of  the  local  contribution  of  the 

E.F.G.  the  value  from  the  graph  in  ( R  a  7  5 )  Fig.  36  was  taken. 

1 7  ? 

This  gives  (1  -  R)eq-|Qc  =  +5.8  x  10  V/cm  .  So  one 

1 7  2 

estimates  eq  =  4.1  x  10  V/cm  ,  yielding  with  =  315 
MHZ  | Q |  =  0.5  barn. 

7  0 

The  temperature  dependence  of  Ga  in  Zn  deviates 
certainly  less  from  the  temperature  dependence  of  qlat  than 
the  various  other  ions  in  Zn  (Fig.  37);  this  might  indicate 
that  the  lattice  contribution  is  in  comparison  larger  than  in 
the  other  cases.  This  would  increase  the  magnitude  of  the 
quadrupole  moment. 

The  other  part,  the  calculated  value  of  Q  is  not 
established.  In  appendix  II  the  results  are  quoted  for  the 
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Fig.  36:  Correlation  of  ionic  and  extraionic  field- 
gradients  in  metals.  Most  values  refer  to 
room  temperature.  A  typical  temperature 
dependence  is  shown  in  the  case  of  ZnCd  by 
an  arc  of  data  points.  Underlined  symbols 
refer  to  the  host  metal.  A  and  B  refer 
to  two  inequivalent  substitutional  sites  in 
o)-Zr.  Filled  circles  indicate  data  from 
e2qQ  values  of  known  signs.  For  open  circles, 
e2qQ  signs  are  unknown.  Locations  of  these 
points  are  predicted.  (Ra75) 
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T(K) 


Fig.  37:  The  variation  of  the  quadrupole  frequencies  of  different  ions 
in  Zn  with  the  temperature  and  the  variation  of  the  lattice 
part  of  the  field  gradient  with  the  temperature. 
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calculation  of  the  energy-levels  in  the  Intermediate  Coupling 

7  0 

Model  for  Ga.  The  results  indicate  that  the  results  can 

be  described  with  some  admixing  of  phonon  states.  With  the 

+  70 

known  quadrupole  moment  of  the  vibrational  2  state  in  Ga 
of  0.38  b  ( C  h  7  2 )  and  the  two  particles  outside  the  core,  a 
quadrupole  moment  of  0.5  b  is  not  unbelievable. 

V 1 . 2  The  Anisotropy  Coefficients 

For  a  determination  of  the  value  of  the  anisotropy 
coefficient  a2  only  the  runs  at  560  K  and  340  K  were  taken 
into  account.  The  three  other  runs  are  measurements  with  two 
angles  only,  which  will  easily  introduce  errors  in  the  values 
large  in  comparison  with  the  runs  with  6  angles. 

The  experiments  yield  different  values  of  a^  for  the 
various  runs.  This  effect  is  not  consistent  between  the  188 
keV  and  the  691  keV  y  rays,  which  indicates  that  misalignment 
is  not  the  only  effect.  Various  other  authors  (Be75,  Va73, 
Bu74,  Ha 7  3)  the  same  effect  when  studying  the  temperature 
dependence  or  pressure  dependence  of  quadrupole  interactions 
with  the  TDPAC  technique.  Two  groups  (Be74  and  Ha73)  report  a 
reduction  of  the  a2  coefficient  by  a  factor  of  two  and  more. 
All  others  show  differences  up  to  25%. 

7  0 

Other  y  rays  of  Ga  showed  agreement  between  the  two 
runs,  and  with  the  anisotropy  coefficients  as  published  by 
( Do7  2  )  and  (Na73). 
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The  theoretical  expression  for  the  time  integrated 
attenuation  coefficient  is  (Fr65) 


Gk(  °°)  EN  skN  1  +  ( Nto  t)2 

q 


(VI. 2) 


where  t  is  the  lifetime. 

The  values  for  this  expression  for  halflife  of  23  ns 
and  spin  4  are  plotted  in  Fig.  38.,  as  a  function  of  w  . 

q 

For  the  560  K  run  the  value  of  the  time  integrated  attenuation 
coefficient  is  0.66  and  the  value  of  w  is  330  MH  .  For  the 

q  z 

340  K  run  we  have  (  00  )  =  0.48  with  =  319  MHz.  If  one  compares 
these  values  with  the  curve  it  can  be  concluded  that  the 
560  K  run  does  not  fit  the  theory  of  a  static  electric  field 
gradient.  The  340  K  run  does  agree  with  the  theory  and  it  can 
be  concluded  that  the  a2(0)  value  is  0.42  ±  0.02  and  a4(0)  = 

0.14  ±  0.05.  For  the  691  keV  y  ray  there  is  no  possibility 
for  such  a  check  for  the  peak  is  obscured  by  the  692  KeV  peak 

7  2 

of  Ge . 

The  three  other  runs,  which  yield  a  mixture  of  a^ 
and  a^,  have  consistent  values  for  the  coefficients  of  both 
transitions,  but  a  temperature  dependence  is  not  excluded  with 
the  present  error  bars. 

If  a  temperature  dependence  is  present  the  results 
might  be  explained  with  the  occurrence  of  a  fraction,  that 
does  not  take  substitutional  lattice  sites,  but  experiences 
a  time-dependent  quadrupole  interaction.  This  fraction  will 
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Fig.  38:  The  time  integrated  attenuation  coefficient  as 

function  of  for  a  level  with  spin  4  and  t^^  = 
23  ns,  for  various  values  of  n  (0.0,  0.2,  0.4, 
0.6,  0.8,  1.0,  with  0.0  the  highest  values  of 
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FREQ  (MHz) 
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increase  with  increasing  temperature.  The  time  differential 
pattern  of  the  anisotropy  can  be  described  with 

G2(t)  =  IN  skN  e  'N6wt  (VI. 3) 

(  where  Sw  is  a  large  frequency  spread  as  in  formula 

(IV.  2)  (Be7  4a,  Bo7  2,  Va73).  This  might  also  explain  the 
results  of  other  results. 

The  constant  c,  which  appears  in  the  fitting  will  in 
these  data  absorb  the  effect  of  the  non-substi tu tiona 1  fraction 
on  the  time  dependent  results.  For  the  340  K  run  it  will 
probably  not  contain  any  effects  of  normalization  errors, 
for  these  errors  will  have  the  same  effect  on  the  time 
integrated  values  as  on  the  time  differential  results  and  the 
present  values  fit  the  theory  for  a  static  electric  field 
grad i ent . 

V  1 . 3  General  Conclusions 

Despite  experimental  difficulties  such  as  the  short 
lifetime  and  a  relatively  small  cross-section,  these  experiments 
yield  the  frequency  and  the  unperturbed  anisotropy  coefficients, 
although  not  the  spin  and  the  n  values. 

The  a2  value  of  0.42  ±  0.02  is  in  fair  agreement  with 
the  value  0.38  predicted  by  the  CN-model,  (Hu75).  It  is  also 
in  agreement  with  the  measured  anisotropy  of  the  411  keV 
( Na 73 ) .  The  a4  value  of  0.14  ±  0.05  shows  the  same  agreement 


as  the  a2  value. 
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The  w  value  of  315  MHz  is  discussed  in  section  VI. 1 
and  indicate  a  quadrupole  moment  of  0.5  b,  when  it  is  assumed 
that  the  universal  correlation  of  the  components  of  the 
field  gradient  as  suggested  in  Ra75  is  true. 

The  TDPAD  technique  has  the  capability  of  measuring 
the  quadrupole  interaction,  but  has  the  drawback  that  it  can 
only  be  performed  with  a  pulsed  beam.  This  gives  us  no 
possibility  to  place  the  nucleus  in  a  particular  place  of 
a  lattice  TDPAC,  which  hasn't  this  drawback,  has  the  dis¬ 
advantage  that  the  level  of  interest  has  to  be  fed  by  a  long- 
lived  parent  (t-j^  >  ^  day).  If  these  techniques  are  performed 
in  combination  wi th  Mos sbauer  and  NQR  experiments,  both  of 
which  have  their  own  drawbacks,  one  is  able  to  measure  field 
gradients  and  quadrupole  moments  of  most  nuclei  in  various 
surroundings.  But  the  difficulties  in  the  measurements  are 
clearly  indicated  by  the  fact,  that  the  theory  was  well 
developed  in  1962,  and  only  in  recent  years  have  the 
experimental  results  been  available  in  reasonable  number. 
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APPENDIX  I 


TARGET  HEATING 

When  a  beam  deposit  its  energy  0  on  the  target, 
radius  b  and  thickness  d,  in  a  circle  with  radius  a  (a  <  b) 
in  the  centre,  the  centre  will  have  a  higher  temperature 
than  the  edge.  The  edge  is  assumed  to  absorb  the  heat  with 
out  a  temperature  rise.  The  temperature  as  function  of  the 
radius  is: 

T(°>  -  T‘ r =  dnr  {ir)2  r  1  a 

T<r>  -  T<a>  -  2^n c  ln  i  a  <  r  <  b 

where  K  is  the  constant  of  thermal  conductivity. 

Define  P  =  0/d,  the  energy  loss  of  the  beam  per  cm  material 
In  this  way 

T  ( b )  -  T  ( 0)  -  ^  U  +  2  In  jf  }. 

Substitute  the  appropriate  values  for  3  MeV  protons  on  a 
zinc  target. 

p  =  70  -  x  7200  mg/cm2  =  504  W  cm"1  pA"1 

mg/cm 2 

K  =  1  . 2  W  /  (cm  K) 
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'  '  U)T  -  (t)T' 
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for  b/a 


3.0 

AT  =  106.9 

2.7 

99.8 

2.5 

94.7 

2.0 

79.7 

1  .8 

72.7 

The  sublimation  rate  of  zinc,  which  was  used  as  a 
check  on  the  approximate  temperature  can  easily  be  calculated 
Every  material  has  a  vapor  pressure,  which  will  be  pumped 
away  in  the  high  vacuum  beam  line.  This  pressure  is  given  by 


1 

3 


P  v 


where  p  is  the  density  of  the  gas  and  v  the  velocity  of  the 
gas  molecules.  Per  unit  time  the  average  distance  travelled 
by  the  molecules  is  v  , .  The  mass  escaping  per  unit  time  and 
per  unit  area  is 


A  m  -  3  P0  vav  /  v 


8m 


7rkT 


{ For  the  formulae  for  vav  and  V7'  see  ref.  1). 

With  the  vapor  pressure  quoted  in  ref.  2  the  sublimation  rate 

70  . 
of  Zn  can  be  calculated  and  the  values  are  tabulated  in 


Table  4. 
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TABLE  4 

SUBLIMATION  RATE  OF  ZINC  AS  FUNCTION 
OF  THE  TEMPERATURE 


T  (K) 

P0  (mm  Hg) 

2 

A  m  (yg/cm  s) 

421 

10"10 

9.7 

X 

10 

449 

io"9 

9.4 

X 

10 

481 

10-8 

9.1 

X 

10 

517 

Hf7 

8.8 

X 

10 

559 

10"6 

8.4 

X 

10 

611 

io"5 

8.1 

X 

10 

ref.  1  D.  Ha  1 1 i day ,  R.  Resnik,  Physics  (J.  Wiley  and  Sons, 
N.Y.,  1966) 

ref.  2  American  Institute  of  Physics  Handbook,  ed .  D.  E. 
Gray,  1972,  p.  4  -  300. 


APPENDIX  II 


THE  LEVEL  SCHEME  OF  70Ga  IN  THE 
INTERMEDIATE  COUPLING  MODEL 

The  Intermediate  Coupling  Model  couples  one  or  two 
particles  in  sheil  model  states  outside  an  even-even  core  to 
the  vibrations  of  the  core.  The  coupling  is  charac ter i zed 
by  the  energy  of  the  phonons,  the  number  of  phonons  and  the 

coupling  constant  £  (Ch67). 

7  0  6  8 
For  Ga  the  even-even  core  is  Zn.  An  estimate  of 

the  phonon  energy  of  this  nucleus  is  1  MeV.  The  particle 

levels  can  be  taken  as  the  single-particle  levels  of  Ni  and 

yield  for  the  energies 

E (  |2p  3/ 2> )  =  0.0  MeV 
E (  | If  5/2>)  =0.78  MeV 

E (  | 2p  1 / 2> )  =1.08  MeV 

E (  | lg  9/ 2> )  =  2.00  MeV 

5  7 

except  for  the  g  9/2  state,  which  is  not  observed  in  Ni.  The 
quoted  value  is  the  estimate  of  Dohan  (Do72). 

When  the  calculations  are  done  with  one  phonon  and 
3  orbitals  available  for  the  proton  and  2  orbitals  for  the 
neutron,  the  results  predict  a  reasonable  part  of  the  observed 
level  scheme.  The  other  levels  might  have  an  admixture  of 
two-phonon  states. 
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The  4  state  at  the  calculated  energy  of  795  keV  above  the 
lowest  state  is  built  up  as: 


0. 

977  l0;  p3/2 

g9/2; 

-0. 

030  |( 

);  p1/2 

99/2; 

-0. 

084  |1 

’  p3/2 

99/2; 

-0. 

079  |1 

;  p3/2 

99/2’ 

+0. 

024  |1 

•  -F 

’  5/2 

99/2; 

-0. 

020  |1 

’  f5/2 

9  9/ 2’ 

+0. 

042  |1 

;  Pl/2 

9  9/ 2’ 

Mi 

f 

T  5/  2 

9  9/ 2  ; 

4> 

i  n 

the  f 

5/2  orbit> 

5/2  y 9/ 2 


3/2  y9/2! 
*3/2  S9/2; 
5/2  g9/2; 
5/2  99/2; 
5/2  99/2; 
’l/2  99/2; 


and 


both  particles  coupled  to  spin  5. 
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Exp. 

E  ( keV)  Jv 


996  — - -  2  + 

901  — — - - 4  + 

879  ———  4- 

691  - - -  2‘ 

651  - - -  1  + 

508  - -  2  + 


Ca!c. 


E  (keV) 

j77 

1  1  10 

.  3“ 

1094  “ 

“  1  + 

998  - — 

~~  0  + 

925 

4  + 

912 

2  + 

475  -  2  + 


130  -  1  + 


0  -  1  + 

Fig.  39:  Comparison  between  the  experimental  and  calculated 

level  scheme  for  ^Ga .  The  calculation  is  done  with 
one  phonon  in  the  Intermediate  Coupling  Model. 
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PROGRAM  TDPAD 


Nuclear  Research  Centre 
University  of  Alberta 
Edmonton ,  Alberta 

INTERNAL  REPORT  #77 

C.W.  Luursema 3  and  D.A.  Hutcheon 


Introduction 

This  program  is  meant  for  the  analysis  of  a  time  dependent  y  ray 
angular  distribution  due  to  static  electric  Quadrupole  Interaction  in 
randomly  oriented  microcrystals  with  the  plane  of  detection  containing 
the  beam.  It  starts  with  time-delay  spectra  of  the  transition  of  interest 
and  of  background  data  for  two  or  more  angles  and  gives  the  least  x2 
fit  parameters  to  the  theoretical  formula. 

It  is  developed  for  the  SDS,  940  for  the  Nuclear  Research  Centre, 
makes  use  of  all  32K  core  and  works  only  with  Rad  Monarch  system. 

TDPAD  is  a  Fortran  Linking  program  consisting  out  of  three  links. 

The  three  links  'Data',  fitting  the  raw  data  to  the  anisotropy  coefficients 
at  various  delay  times,  'Display',  fitting  the  coefficients  to  the  theoretical 
expression  and  central  processor  and  'Sun'  calculating  the  coefficients, 
are  reasonably  independent  so  the  program  doesn't  switch  very  often  during 

running. 

It  is  assumed  that  the  raw  data,  i.e.  the  yield  as  function  of  delay, 
are  stored  in  maximal  1024  channels  with  the  time  running  backwards,  i.e. 
increasing  channel  number  corresponds  to  decreasing  delay. 
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Theory 
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Time  dependent  Perturbed  Angular  Distribution  is  given  by  (all  refs.) 


W(e>)  -  ^(a^  (t)  +  c)  PjJcoso) 
where  c  is  a  time  independent  constant. 

Gk  (t)  =  Eqmm"  21c?T  (I  -■  I  m'|k  q)  * 
snn,  (I  -m"  I  m" 1 | k  q)  e'i(En 


x 


u  u*  „  u  u*  ,  , 

mn  m  n  m  n  m  n 1 

where  U  is  the  transformation  matrix  which  diagonalizes  the  Hamiltonian 
(UHU-1  =  E), 

Since  (E  -E  ,)  =  -(E  ,  -  E  )  we  end  up  with  only  cosine  terms 
nn  n  n  r 

G|(,(t)  =  ZN  skN  cos  Nut 


-  (E„  -  E  , ) 


where  o>  =  e2qQ/h  and  N  =  v  n  n 

ezqQ 

The  Hamiltonian  is  given  by  the  elements 

eQ 


H 


m 


(3m2  -  I  (1+1 ))  V 


zz 


H 


eQ 


m,m±2  =  n  4I(2I-lT  X 

{(  I  +  m-1)  (I  +  m)  (I  ±  m  +  1)  (I  ±  m  +  2)} 

with  V  =  - —  V,  evaluated  at  the  nuclear  site 
zz  az2 

V  -  V 

— and  0  £  n  l  1 
zz 


and  n  = 


' 
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By  taking  into  account  a  spread  of  frequencies  around  the  central 
one  and  a  finite  resolution  time  of  the  detecting  system,  the  formula 

Gk(l)  =  E  SkN  cos  N  „t  e'5“Nt  e'1/2(“NT)2 

6w  =  absolute  frequency  spread 
x  =  resolution  time  of  system 
The  fitting  is  done  to 

al<(t)  =  ( t )  ak(0)  +  c, 

where  c  is  a  time-independent  constant. 

Procedure 

The  program  is  split  up  in  three  parts 

(a)  'Data1.  This  part  reads  in  the  data  from  magnetic  tape  unit  number  1 
and  calculates  the  a^  coefficients  as  a  function  of  time.  The  input 
data  consist  of  the  time  spectrum  for  the  region  of  interest  and  the 
appropriate  background  spectrum,  both  should  be  in  the  same  run  on 
the  tape.  It  is  possible  that  there  are  two  angles  in  one  run. 

After  reading  in  data  for  a  least  two  angles,  it  orders  the  data 
into  groups  (maximum  of  512)  and  fits  each  group  to 
W(o)=EkakPk  (cose)  for  K=0,2,4. 

If  desired  a  list  of  calculated  values  and  their  errors  is  printed 
out  (time  runs  backwards  on  this  list). 

The  search  for  next  runs  is  according  to  run  number,  i.e.  a  higher 
run  number  causes  a  search  forward  of  the  tape  and  a  lower  search 
backwards.  By  reading  in  a  faked  high  or  low  run  number  with  a  blank 
card  behind  it,  it  is  possible  to  direct  the  search  freely  without 
changing  the  data. 
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(b)  'Display'  The  central  processing  and  fitting  unit.  This  part  asks 
for  the  other  information  necessary  for  fitting  to  the  formula. 

(1)  The  out  put  from  'Data',  i.e.  a^(0)  and  the  error  in  a^(0). 

After  calculation  of  ak(0)  and  x2  it  can  be  directed  to  all 
parts  of  the  program  (see  listing  in  Table  I).  The  program 
types  after  each  step  the  number  of  channels  used  for  the  fit, 
the  a^(0)  value,  its  the  error  in  ak(0)  and  the  x2- 

(c)  'Sun'  This  part  calculates  the  S  ^  coefficients  and  the  values 
of  N  for  formula  (1)  as  function  of  spin,  n  and  k.  It  calculates 
the  matrix  and  diagonalizes  it  with  the  method  of  Jacobi  then  it 
calculates  for  each  N,  making  use  of  the  fact  that  for  each  n 
only  half  of  the  magnetic  substates  contribute. 

Table  I  List  of  options  in  'Display'  on  question  'O' 

0.  Types  out  all  possibilities. 

1.  Goes  to  'Data'  and  read  in  new  set  of  data. 

2.  Goes  to  'Sun'  and  asks  for  new  values  of  k,  spin  and  n  and  recalculates 
SkN  and  N. 

3.  Print  out  all  relevant  parameters,  e.g.  S ^  coefficients,  N-values, 
ak(0),  error  in  ak(0),  x2>  spin,  n,  k  first  and  last  channel  of 
the  region  of  fitting,  time-zero,  frequency,  spread  in  frequency 
and  fudge. 

4.  Asks  for  new  value  of  time-zero. 

5.  Asks  for  new  region  of  fitting. 

6.  Asks  for  a  new  frequency,  constant  term  and  relative  spread. 

7.  Goes  to  show-routine,  which  shows  experimental  data  with  upper  and 
lower  limit  of  the  error  bar,  the  calculated  curve  and  the  zero  line 


n  log'  no  Jwfr  tost  etti  to  szh  g/nift*  t  .1  does  *<ot  .jt  *etormfh3 
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if  the  upper  and  lower  limits  (amax  and  amin)  contain  zero.  To  leave 
this  routine  set  breakpoint  4. 

8.  Goes  to  the  plotting  routine,  if  and  only  if  the  program  has  displayed 
the  data  after  the  last  calculation,  otherwise  it  does  the  same  as 

7. 

9.  Goes  to  the  plotting  routine,  taking  old  values  for  the  plotting 
parameters  except  for  YOFF,  unless  they  are  zero,  in  which  case 
it  goes  to  7. 

Input  for  Data 

Card  I:  NGRP,  NCPG,  NLEG,  ID1R,  TITLE  (FORMAT  414,  10A4). 

NGRP  number  of  groups  less  than  513. 

NCPG  number  of  channels  per  group  NGRP  *  NCPG  <_  1024. 

NLEG  number  of  Legendre  polynomials  (maximum  3). 

IDIR  direction  first  search  (0  forward,  1  backward). 

Title  (less  than  40  characters). 

Card  II:  Date,  Month,  Year  Run,  Angle 
Format  413,  F6.0 . 

Year  =  0  calculations  follows. 

Card  III:  IBASE,  ZNORM 
FORMAT  15,  F7.3 

IBASE:  First  Channel  of  spectrum  to  be  read. 

ZNORM:  Normalization  constant  with  which  the  data  are  multiplied 
before  processing.  ZNORNkO  Background 

ZN0RM>0  Real  data 


ZNORM  =  0  go  to  next  angle. 
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So  proper  order  is: 
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I,  II,  III  (background),  III(real  data),  Ill(blank),  II,  Ill(background) , 
III (real  data) , I Il(blank) ,  Il(blank) 

Teletype  input 

PRINT  YES ( 1 )  OR  NO  (0)  (II) 

1:  PRINTS  out  ALL  a^  for  each  group  normalized  to  aQ=l  and 

not  normalized  with  group  number,  error  in  a^  and  x2  of 
fit  for  each  group. 

0:  Suppresses  this  output,  and  this  information  is  lost,  except 

for  the  nomalized  values,  which  are  used  by  'DISPLAY'  and 
may  be  plotted. 

ENTER  1  TO  RESTART,  0  TO  EXIT  (II) 

1:  Begins  again  at  the  beginning  of  the  program  (reread  all  data) 

0:  Stops 

Resolution  Time  (FWHM) 

Floating  point  input  of  the  resolution  time  of  the  electronics 
and  detector  expressed  in  number  of  groups. 

FREQ  IN  RAD/GR,  C,  REL  SPREAD 
Format  3F10.0 

Asks  for  the  essential  parameters  in  the  formulae,  the 
frequency  being  wq  and  relative  spread  being  6^ 

TZ(F  10.0)  M 

Input  of  time  zero,  as  group  number  according  to  the  output 

1 i sting  of  'Data ' . 

AMIN,  AMAX 

The  lower  and  upper  limit  of  the  values  for  the  show-routine 
and  for  the  plot  routine.  Values  outside  this  region  are  set 
to  the  upper  or  lower  limit. 


. 
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GROUPS/INCH,  YSC,  YOFF  (Floating  point) 

GROUPS/INCH  X-scale 

YSC,  number  of  inches  over  which  the  y-scale  extends. 

YOFF,  vertical  offset  of  AMIN  in  inches,  zero  being  point  of 
pen  location. 

EXP(l)  CAL (2 ) ,  AX IS ( 3 ) ,  SLEW(4),  D0NE(5) 

Format  II 

1.  Plots  experimental  error  bars. 

2.  Plots  the  current  theoretical  curve. 

3.  Plots  the  x-and  y-axis,  from  time-zero  to  the  last  channel  to 
plot  and  the  y-axis  from  amin  to  amax,  and  x-axis  at  a^=0  and 
the  y-axis  on  time  =0. 

4.  Slews  the  paper  up  till  the  second  integer  inch  after  the  plot 
and  cannot  be  returned. 

5.  Returns  to  the  rest  of  the  program. 
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Plot  same  as  above 

Other  questions  are  self-explanatory. 

All  input  parameters  are  floating  point  unless  otherwise  indicated. 
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